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(short  pulse)  radar.  Theoretical  details  pertaining  to  this  approach  are 
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representing  differently  shaped  regions  on  the  target  surface.  The 
dimensionality  of  the  feature  vectors  thus  constructed  is  very  low  (less 
than  ten).  Moreover,  when  used  in  the  classification  of  targets  by  a 
nearest  neighbor  classification  strategy,  such  feature  vectors  permit 
accurate  discrimination  between  targets  that  do  not  differ  much  in  shape; 
and  also  they  are  in  a  large  measure  insensitive  to  noise. 

The  above  results  were  corroborated  by  computer  simulations 
performed  on  the  data  base  created  by  the  coherent  X-band  short  pulse 
(0.5  nanosecond)  radar  at  the  Fort  Worth  operation  Radar  Range  of  General 
Dynamics  Convair  Aerospace  Division.  The  three  objects  on  which 
extensive  simulations  were  carried  out  were:  in  AGENA  space  vehicle  with 
rectangular  cross-sectional  first  stage,  an  AGENA  space  vehicle  with 
cylindrical  cross-sectional  first  stage  and  an  AGENA  space  vehicle 
payload.  The  results  of  these  simulations,  which  were  rather  sucessful, 
are  discussed  in  detail. 


[1]  R.J.P.de  Figueiredo  and  C.L.  Ha,  "A  Four ier-Prony  Tauberian 
Approach  to  the  Analysis  of  a  Mixture  of  Delayed  Signals",  Rice  University 
Technical  Report  EE-7907,  Houston,  Texas,  September  1979. 
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1.  INTRODUCTION 


In  a  wide  bandwidth  (short  pulse)  radar  system*,  the  return  signal 
consists  of  a  sequence  of  short  pulses  emanating  from  different 
scattering  centers  on  the  target  [2,  3,  4,  5].  These  centers  represent 
differently  shaped  regions  on  the  target  surface.  Figure  1  illustrates 


Range  Gate 


Fig.  1  Typical  Short-Pulse  Signature 


the  X-band  radar  return  signature  from  1/8-scale  model  of  an  AGENA 
vehicle  taken  from  reference  [3].  The  first  pulse  is  the  return  from  a 
fixed  reference  scatterer,  while  the  remaining  pulses  occur  at  the  points 
of  discontinuities  on  the  vehicle  surface  as  shown. 


*Wide  bandwidth  radar  systems  are  of  recent  vintage  and  possess  the  very 
high  resolution  needed  to  identify  complex  objects.  Examples  of  such 
systems  are  the  SPTF  facility  at  Floyd,  New  York  and  the  ARPA-Lincoln 
ALCOR  facility  at  Kwajalean. 
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In  what  follows,  we  use  a  Tauberian  representation  for  the  return 
signal  y,  that  is 


M 

y(t)  =  T  a .  -  O,  (D 

i  =  l 

where  x  is  the  reference  pulse,  shaped  as  the  incident  pulse,  and  aj 
and  Tj  are  the  reflection  coefficient  and  delay  associated  with  the  ith 
scatterer.  In  the  more  general  case  in  which  the  scatterers  are 
distributed  over  the  larger  surface,  we  have  a  representation  for  y  in 
the  form 


where 


y(t)  =J*  h(t  -  t')  x(t')  dt',  (2) 

M 

h(t)  =  Z  a.  »(t  ■  t.  -  t  )  (3) 

I  I  U 

1=1  '  ' 


is  the  inpulse  reponse  of  the  target;  and  is  the  delay  corresponding 
to  the  target  range. 

Let  the  data  consist  of  a  set  of  samples  x  £  x(tn  +  n  At)  and 

n  u 

y  =  y(t  -  tq  +  n  At),  n  =  0,  1,  ...,  N-l,  of  x  and  y,  where  the 

constants  t^,  T^f  and  At  are  conveniently  chosen.  In  the  present  paper, 

we  propose  the  use  of  a  frequency  domain  Prony  approach,  discussed  in 

detail  elsewhere  [1],  to  determine  the  parameter  vectors  _t  =(Ti,  . . . ,  TM  ) 

and  £  =  (a  ,  ...,  a  ),  from  the  above  data.  One  can  then  use  an 
l  M 

appropriate  subset  of  the  variables  a]}  V  a2’  V  • • • >  V  TM  »  and  of 
the  parameters  in  h  5,  as  the  components  of  a  feature  vector  w  in  a 
nearest  neighbor  algorithm  for  target  classification.  In  the  present 
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paper,  we  limit  our  attention  to  the  representation  (1),  the 
generalization  to  the  case  of  equations  (2)  and  (3)  being  clear  from  the 

developments  in  [1].  } 

| 

In  section  2,  we  briefly  describe  the  frequency  domain  Prony 
approach  used  here.  In  sections  3  and  4,  we  present  the  results  of  some 
of  the  simulations  performed  on  the  signature  data  obtained  through  use 
of  the  coherent,  X-band  short-pulse  (0.5  nanosecond)  radar  at  the  Fort 

j; 

Worth  operation  Radar  Range  of  General  Dynamics  Convair  Aerospace 

| 

Division.  The  data  was  provided  to  us  on  tape  by  the  Rome  Air 

J  f 

Development  Center  (RADC)  (Air  Force  Systems  Command).  Specifically  we 
are  concerned,  in  section  3,  with  the  properties  of  the  Tauberian 

representation  (1),  and,  in  section  4,  with  the  application  of  the 

,  | 

nearest  neighbor  algorithm  for  target  classification  using  a  few  of  the  i 

| 

components  of  <3  and  jr  as  features.  Using  in  some  cases  not  more  than  8 
features  and  not  less  than  four  features,  accurate  discrimination  among 
three  space  vehicles  was  achieved. 

2.  THE  FREQUENCY  DOMAIN  PRONY  ALGORITHM 

! 

j: 

The  application  of  the  frequency  domain  Prony  algorithm  for  analysis  ! 

of  the  return  signature  for  short  pulse  radar  was  first  considered  by  H. 

Webb  [6]  and  E.  M.  Kennaugh  [7].  Frequency  domain  Prony  methods  for  both 

noise-free  and  noisy  data  were  investigated  in  detail  by  the  authors  in 

[1J.  j 

To  apply  the  Prony  technique  to  equation  (1),  we  simply  take  the 
Fourier  transform  of  (1)  to  obtain 
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h(u>)  =-l(y)a 
Tour 


M  "j  T(  (1) 
E  aie 
i  =  l 


(4) 


From  the  smaples  xQ,  x^,  and  yQ,  y^  ,  one  calculates  by 

DFT  a  set  of  samples  of  H  at  equidistant  frequency  points  uu  =  uj  +  kAuu  , 

k  0 

k  =  0,  . ..,  K  -  1,  where  2M  £  K  £  N. 

Let  us  use  the  notation 

Hk  =  H(ojk).  (5) 

We  use  the  Prony  algorithm  to  find  the  2M  parameters  in  the  vectors 
r  =  (  T  t  . . . ,  t  )  and  a  =  (a  ,  . . . ,  a  )  by  fitting  the  complex  numbers 

—  IN  IN 

Hk,  calculated  from  the  data,  to  the  expression  (4),  with  equality  or  in 
the  least  squares  sense,  i.e., 


M  "J  Ti 

E  a.  e  =  H  ,  k  =  0,  1,  K-l. 

i=»l  1  k 


(6) 


A  modification  of  the  conventional  Prony  method  to  conveniently  take 

care  of  the  fact  that  H  are  complex- valued  (rather  than  real  valued  as 

k 

in  the  conventional  case)  is  presented  in  [1],  We  have  also  shown  in  [1] 
how  to  apply  our  procedure  when  the  smaples  are  noisy  and  also  how  to 
take  care  of  the  general  case  of  the  expressions  (2)  and  (3)  using  Prony 
techniques.  In  the  present  paper  we  do  not  discuss  these  developments 
and  simply  assuming  that  the  return  signal  has  the  form  (1),  and  given 
the  samples  of  x  and  y,  we  determine  and  by  the  procedure  mentioned. 

We  now  proceed  to  explain  our  approach  to  radar  signature 
classification  in  the  context  of  simulation  results  performed  both  on 
noiseless  and  noisy  data. 
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3.  TAUBERI AN  RECONSTRUCTION  OF  THE  RETURN  SIGNAL 


The  actual  radar  return  signal  y  consists  of  samples  of  a  frequency 
carrier  in  the  X-band,  single-sideband  modulated  by  the  scatterer 
interaction.  Hence  J  has  the  form 

y(n)  =  ct( t Q—  t0+  n  At)cosw0(t0-  to+  n  At) 

+S(to-T0+  n  At)sinwQ(to-To+  n  At) 

n  =  0,  1,  ...t  N-l ,  (7) 

where  u>Q  denotes  the  carrier  frequency,  and  <*(.)  and  6(.)  are  functions 
representing  the  modulation.  By  passing  y  through  a  synchronous  detector 
system  such  as  the  one  shown  in  Figure  2,  the  envelope  and  phase 
information  in  terms  of  the  variables  yp  and  0n  are  defined  by  (8a)  and 
(8b  j; 


2 


complex  sampled 
data  signal 


t 


n 


a 

n 


P 


n 


=  -  tq  +  "At 


Fig.  2  Synchronous  Detection  System 
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y  =  y(t  “T.+  n  At)  =  faxt  -t  +  n  At)+82  (t  -t  +  n  At),  (8a) 

n  0  0  -I00  00 

9  =  ©( -t  +  n  At)  =  tan  [8(t  -  t  +  n  At)/fct(t  -t  +  n  At)], 

n  u  0  0  0  0  0 

n  =  0,  N-l  (8b) 

The  phase  sequence  (  0  ]  was  not  used  in  our  classification 

n 

algorithm.  For  each  space  vehicle  model  investigated,  the  envelope 

sequence  fy  ]  for  a  given  yawing  angle  was  obtained  from  tapes  supplied 
n 

to  us  by  RADC.  The  reference  pulse  (the  first  pulse  in  Figure  1)  was 
used  as  the  waveform  of  x. 

Our  first  set  of  simulations,  described  in  the  present  section,  were 
designed  to  test  the  representation  (1)  by  reconstructing  the  actual 
waveform  y  as  the  waveform  9  described  by 


9(n) 


M 

E  a .  x(n  -  t.  ) 
i  =  l  '  1 


(9) 


where  the  constants  a.  and  t.  were  calculated  by  fitting  the  data  f  yn  :  n 

=  0,  N-l]  and  (xn:  n  =  0,  N-l]  to  the  expression  (1)  by  our 

frequency  domain  Prony  technique,  i.e.  we  calculated  H.  by 

k 


N-l 

E 

n=0 

N-l 

E 


Y 

n 


X 


-jw 

e 


k 


n 


k=  0,  1,  ...,  K -1  (10) 


and  then  obtained  jt  and  £  by  satisfying  (6)  using  the  Prony  method.  The 
number  M  of  delays  was  set  equal  to  an  integer  slightly  greater  than  the 
number  of  peaks  of  high  amplitude. 

The  data  used  is  described  in  the  reports  by  Cisco,  Johnston,  and 
Gruver  [2],  Cisco,  Coble,  and  Gruver  [3], 

In  our  simulations  we  selected  as  targets  the  following  three 
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objects  at  yawing  angles  indicated: 

Model  59:  AGENA  space  vehicle  with  rectangular  cross-sectional 
first  stage,  yawing  angle  6=  0.2°; 

Model  66:  AGENA  space  vehicle  with  cylindrical  cross-sectional 
first  stage,  yawing  angle  8=  0.2°; 

Model  60:  AGENA  space  vehicle  payload,  yawing  angle  0  =31.2°. 

For  each  case  the  demodulated  radar  return  signal  recorded  in  the 
tape  is  exhibited  in  Figures  3,  4,  and  5  respectively.  Figures  6,  7,  and 
8  show  the  locations  of  these  signatures  in  3-dimensional  plots  of 
intensity  versus  time  and  yawing  angle  (obtained  from  [2]). 

We  now  comment  on  the  reconstruction  of  the  return  signature  for 
each  of  the  above  models  separately. 

*  Model  59: 

We  first  removed  the  bias  in  the  signature  for  this  model  (appearing 
in  Figure  3) and  selected  for  display  in  Figure  9  the  useful  sections  of 
the  signature,  namely  the  reference  signal  x  and  the  return  y,  side  by 
side. 

We  assumed,  first,  the  order  (number  of  delays)  M  of  the  system  to 
be  4,  and  applied  the  frequency  domain  Prony  method  to  this  data. 

The  signal  y  reconstructed  according  to  (9)  using  the  parameters  t 
and  a  thus  calculted  is  shown  by  the  full  curve  of  Figure  10  (the  true 
signature  y  is  represented  by  the  dotted  curve  in  this  figure).  The 
negative  peak  results  from  the  physically  unmeaningful  fact  that  the 
parameter  83  calculated  above  is  negative.  We  concluded  therefore  that 
the  third  delay  was  redundant  and  reduced  M  to  3,  and  applying  the  same 
procedure,  obtained  the  values  for  t  and  a  tabulated  in  Table  1  in  the 
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Fig.  3  The  Reference  Signal  x  and  the  Return  y 
for  Model  59  at  Yawing  Angle  0.2  degree 
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Fig.  4  The  Reference  Signal  x  and  the  Return  y 
for  Model  66  at  Yawing  Angle  0.2  degree 
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Fig.  5  The  Reference  Signal  x  and  the  Return  y 

for  Model  60  at  Yawing  Angle  31.2  degrees 
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Response  of  Aerospace  Model  66 


Fig.  8  Response  of  Aerospace  Model  60 
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Fig.  10  Comparison  of  the  Actual  Observed 
Signature  with  the  Reconstructed 
Signature  for  Model  59  at  Yawing 
Angle  0.2  Degree 
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row  labeled  "Model  59:  clean  data". 

The  waveform  reconstructed  by  means  of  these  parameter  vectors  is 
shown  by  the  full  curve  of  Figure  11.  We  note  that  it  agrees  very 
closely  with  the  actual  signature  shown  by  the  dotted  curve. 

Model  66  and  60: 

The  delay  and  coefficient  vectors  _t  and  a  were  calculated  in  the 
same  way  as  for  the  model  59,  they  are  listed  in  the  rows  of  Table  1 
labeled  as  clean  data  for  these  two  models.  The  corresponding 
reconstruction  y  of  the  signal  y  is  shown  respectively  in  Figures  12  and 
13. 


Table  1 


Clean 

Delays 

11.4765 

24.9251 

47.1603 

Model 

Data 

Amps 

0.3699 

0.4304 

0.3935 

59  a£ 

Noisy 

Delays 

11.3001 

25.3624 

47.1058 

fl=0.2° 

Data 

Amp 

0.3794 

0.5373 

0.3673 

Model 

Clean 

Delays 

5.0879 

24.4645 

'37.7121 

55.9327 

Data 

Amps 

0.1661 

1.7065 

0.3537 

1.2050 

60  at 
0=31.2° 

Noisy 

Delays 

4.4692 

24.6126 

38.0665 

56.5352 

Data 

Amps 

0.0177 

1 . 8748 

0.3418 

1.2164 

Clean 

Delays 

17.4324 

23.1708 

Model 

Data 

Amos 

0.5375 

0.9985 

66  at 

Noisy 

Delays 

17.3490 

22.9287 

0=0.2° 

Data 

Amps 

0.6203 

0.9188 

In  order  to  test  the  robustness  of  our  procedure,  we  contaminated 
the  data  with  additive  white  Gaussian  noise  with  zero  mean  and 
signal- to-noise  ratio  equal  to  7.75.  The  results  of  the  application  of 
our  procedure  modified  for  the  noisy  case  according  to  the  developments 
in  [1]  are  shown  in  the  rows  of  Table  1  labeled  as  pertaining  to  noisy 
data.  We  note  that  there  is  very  little  change  in  the  parameter  values 
due  to  noise. 


1AX.  RANGE  *  0.  1006*02 
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Comparison  of  the  Actual  Observed 
Signature  with  the  Reconstructed 
Signature  (  M  =3  )  for  Model  59 
at  Yawing  Angle  0.2  Degree 


Fig..  11 


Fig.  13  Comparison  of  the  Actual  Observed 
Signature  with  the  Reconstructed 
Signature  for  Model  60  at  Yawing 
Angle  31.2  Degrees 
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For  model  59,  we  further  investigated  the  effect  of  noise  on  the 
calculated  parameter  values  for  different  signal-to-noise  ratios.  The 
results  are  summarized  in  Table  2.  Again,  there  is  very  close  agreement 
between  the  corresponding  noise-free  and  noisy  cases.  These  remarkably 

Table  2 


■aii»ifSBET itm 

f«3ETC» 

24.9251 

47.1603 

Noise 

HH 

0.3699 

0.4304 

mmmm 

S/N=7. 75 

mrwvim 

25.3624 

0.3794 

0.5373 

0.3673 

S/N=3. 88 

kNSssBi 

11.0809 

24.7749 

47.5632 

0.3038 

0.3922 

0.3039 

S/N=2. 58 

10.9414 

24.8378 

47.7400 

■ES&SMi 

0. 2947 

0.3922 

0.2801 

S/Nxl.94 

Delays 

10.8262 

24.9027 

47.9388 

wmfuZWmm 

0.2847 

0.3979 

0.2545 

S/N=1.29 

iViTSI.l'i.ljMI 

■ ITOREWi 

25.0273 

48.3992 

Amps 

0.2658 

0.4122 

0.2000 

good  results  are  highlighted  by  the  curves  of  Figure  14  and  15.  In 
Figure  14  the  reconstructed  signature  y  with  parameters  calculated  from 
noisy  data,  for  signal-to-noise  ratio  =  1.29,  is  compared  to  the  actual 
noisy  signature  (i.e.  original  signature  with  noise  added).  Figure  15 
compares  the  same  signatue  9  (reconstructed  from  the  noisy  data  with 
signal-to-noise  ratio  =  1.29)  to  that  reconstructed  from  clean  data. 

These  simulations  clearly  attest  to  the  good  robustness  of  our 


procedure. 


Z  Z 
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4.  A  PATTERN  RECOGNITION  ALGORITHM  FOR  RADAR  RETURN 
SIGNATURE  CLASSIFICATION 


Suppose  there  are  n  standard  models  whose  return  signatures  are 

known.  Let  and  a.*  denote  the  parameter  vectors  and  a  belonging  to 

the  i_th  model.  Then  we  denote  by  M  the  dimension  of  the  longest  vector 
I 

1  i  i  =  •••»  nJ  and  we  change  the  lengths  of  all  the  vectors  t*  and  a' 

to  M  using  the  convention  that  if  the  kth  model  has  t  of  length  M  <  M 

k  ^ 

then  the  first  Mk  entries  of  will  consist  of  the  components  of  t_  for 

this  model  and  the  remaining  entries  will  be  set  equal  to  zero.  The  same 

rule  applies  of  course  to  £k. 

Since  the  relative  differences  between  the  components  of  t1  and  the 
relative  ratios  between  the  (nonzero)  components  of  a  1  remain  invariant 
for  a  given  model  i,  we  form  two  feature  vectors  t1  and  a  '  defined  by 


and 


~  l  i  i  . 

X  '  X|  J  •  •  •  9  T  ^)  , 


—  -  »  3  ^  ) 


I  i  i 

Tk=  Tk-  Ti  *  i  -  \  * 0 

~i  i  . 

a .  =  a i.  /(0.1  max  a  .  ) 

k  j  J 

(where  and  ak  denote  the  parameters  Tk  and  ak  for  the  ith  model 
defined  in  the  preceding  sections). 

A  simulation  was  performed  in  which  the  values  of  the  feature 
vectors  r\  and  a  \  i  =  1,  2,  3,  for  the  three  models  59,  60,  and  66 
obtained  from  clean  data  in  the  preceding  section  (Table  1)  were 
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considered  as  "trained  models". 

Then  four  incoming  model  signatures  were  fabricated  by  adding 
different  amounts  of  zero  mean  white  Gaussian  noise  to  the  signatures  of 
models  59,  60,  and  66.  Specifically  the  four  models  were  created  with 
(1)  model  59  noisy  signature  with  signal-to-noise  S/N  =  7.75;  (2)  same 
as  (1)  except  that  S/N  =  1.29;  (3)  Model  60  noisy  signature  with  S/N  = 
7.75;  and  (4)  model  66  noisy  signature  with  S/N  =  7.75.  The  values  of 
the  feature  vectors  i[  and  la  for  each  of  these  signature  were  calculated 
by  our  method. 

A  nearest  neighbor  algorithm  was  used  to  classify  each  of  the 
incoming  models.  For  this  purpose,  the  Euclidian  distance  between  the 
values  of  the  j  vectors  was  used  as  the  distance  measure.  For  those 
trained  models,  for  which  certain  components  of  the  vectors  J_  were  very 
close,  we  used  the  values  of  the  corresponding  components  of  the  vector]* 
as  some  of  the  feature  variables  (in  addition  to  the  delay  feature 
variables  mentioned  above).  This  nearest  neighbor  algorithm  classified 
correctly  all  the  above  four  incoming  signatures. 


5.  CONCLUSION 

The  preceding  results  show  that  the  signal  processing  of  wide 
bandwidth  radar  returns  by  frequency  domain  Prony  methods  constitute  a 
very  economic  way  of  reducing  the  dimensionality  (extracting  features)  of 
the  signatures  observed,  and  lead  to  classification  algorithms  which  are 
very  accurate  even  in  the  presence  of  considerable  noise.  These 
developments  clearly  extend  to  the  case  of  equations  (2)  and  (3)  as 
discussed  in  [1]. 
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